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Maternal obesity and diabetes 
induces latent metabolic defects and widespread 
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Intrauterine nutrition can program metabolism, creating stable changes in physiology that may have significant health 
consequences. The mechanism underlying these changes is widely assumed to involve epigenetic changes to the 
expression of metabolic genes, but evidence supporting this idea is limited. Here we have performed the first study of 
the epigenomic consequences of exposure to maternal obesity and diabetes. We used a mouse model of natural-onset 
obesity that allows comparison of genetically identical mice whose mothers were either obese and diabetic or lean with a 
normal metabolism. We find that the offspring of obese mothers have a latent metabolic phenotype that is unmasked by 
exposure to a Western-style diet, resulting in glucose intolerance, insulin resistance and hepatic steatosis. The offspring 
show changes in hepatic gene expression and widespread but subtle alterations in cytosine methylation. Contrary to 
expectation, these molecular changes do not point to metabolic pathways but instead reside in broadly developmental 
ontologies. We propose that, rather than being adaptive, these changes may simply produce an inappropriate response 
to suboptimal environments; maladaptive phenotypes may be avoidable if postnatal nutrition is carefully controlled. 



Introduction 

Obesity and its accompanying morbidities, such as type 2 diabe- 
tes, are now pandemic. While there is no doubt that obesity has 
a genetic component, its rapidly increasing prevalence indicates a 
major environmental contribution; the etiology and familial pat- 
terns of obesity and diabetes are thus of great interest. Extensive 
evidence indicates that maternal nutrition and metabolism can 
stably affect the phenotypes of offspring, a phenomenon termed 
fetal programming. 1,2 Rodent models of maternal overnutri- 
tion — diet-induced obesity and diabetes (through high fat and/ 
or sugar diets), genetic manipulation (e.g., ob or db mice) and 
chemical/surgical induction of diabetes 1,3,4 — show a variety of 
metabolic consequences in offspring, including increased body 
weight, adiposity, hypertension, insulin resistance, hyperglyce- 
mia, hyperphagia and leptin resistance (reviewed in refs. 1 and 
3). Of these models, diet-induced maternal obesity and diabetes 
is most likely to reflect the metabolic syndrome seen in human 



populations but cannot distinguish the effect of excess fat con- 
sumption from the effect of maternal obesity and diabetes per 
se. Increasing the fat content of a diet also necessitates a relative 
reduction of other components in the diet. Genetic 5 and pairwise 
high fat-feeding models 6 have suggested that maternal obesity per 
se may cause programming effects that are distinct from those of 
dietary fat intake. 

Since fetal programming appears to involve stable changes in 
phenotype, it has been widely speculated that epigenetic modi- 
fications induced by the intrauterine environment are respon- 
sible. 1,2,7 Epigenetic modifications specify stable states of gene 
expression; they are laid down largely during early development 
and remain relatively fixed over a lifetime (and sometimes into 
succeeding generations). In this view, fetal programming stems 
from induced epigenetic changes in genes that regulate metabolic 
pathways and has an adaptive function. But evidence for such 
targeted epigenetic change is scarce, and studies coupling the 
physiological and epigenomic consequences of exposure to fetal 
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overnutrition, particularly in models of natural-onset obesity 
and diabetes, are lacking. It is at present unclear to what extent 
maternal obesity impacts on offspring phenotype via changes to 
epigenetic states. 

We have previously characterized changes in cytosine meth- 
ylation patterns in isogenic C57BL/6J mice in response to dietary 
intervention. 8 Cytosine methylation is a key epigenetic modifica- 
tion in vertebrates, acting as a focal point for mechanisms that 
suppress gene expression. 9 We found spontaneous and wide- 
spread methylation variability in genes with functions related 
to gene expression and development; similar variation has also 
been observed in humans 10 and methylation at four such meth- 
ylation-variable loci appear to correlate with body mass index. 11 
We found that exposure to dietary methyl donors dramatically 
increases the extent of this epigenetic variability. 8 We therefore 
wondered whether similar epigenomic effects would be manifest 
with exposure to maternal obesity. 

Here we have studied the phenotypic and epigenomic effects 
of maternal obesity in a model of natural-onset obesity and type 
2 diabetes, the A" y mouse. Our system allows comparison of 
genetically identical mice whose mothers either did, or did not, 
acquire diabetes and obesity eating a standard chow diet. We 
find that the offspring of obese mothers have a cryptic pheno- 
type that is unmasked by exposure to a Western-style high fat 
diet: they remain nearly normal on a standard diet, but on a 
high fat diet they develop insulin resistance, glucose intolerance 
and hepatic steatosis. This predisposition is associated with gene 
expression and epigenomic changes affecting a broad variety of 
functions. 

Results 

A model of natural-onset maternal obesity. To observe the phe- 
notypic and epigenomic consequences of natural-onset maternal 
obesity and diabetes, we exploited the obese yellow phenotype 
exhibited byA" y mice. 12,13 TheA" y allele is unusual in that it is sub- 
ject to stochastic epigenetic silencing in an isogenic background, 
so that genetically identical mice can display markedly different 
phenotypes. Mice carrying an active allele have a characteristic 
syndrome of fully yellow pelage, obesity and type 2 diabetes; 
mice carrying a silent allele have agouti fur and a metabolically 
normal phenotype, termed pseudoagouti. 14 The A" y allele is con- 
genic with nonagouti (a), and A"? mice are typically maintained 
a.sA" y /a heterozygotes; in a cross of a.i\A" y la dam with an ala sire, 
50% of offspring will be heterozygous for the A" y allele, and the 
other 50% will be homozygous a/a, on an isogenic C57BL/6J 
background. 

In this experiment, we compared the ala offspring of obese yel- 
low diabetic A" y /a dams with genetically identical ala offspring of 
lean pseudoagouti A° y la dams (Fig. 1A). In this way, phenotypic 
and epigenomic parameters could be assessed without confound- 
ing genetic variation. Furthermore, the induction of obesity and 
diabetes in yellow A" y mice does not involve any external manipu- 
lation: the mice are hyperphagic and by 3 mo of age they are 
obese, demonstrably hyperglycemic, hyperinsulinemic and insu- 
lin resistant; thus yellow coat color is inseparable from type 2 



diabetes in adult A" y mice. 13,15 " 17 We chose the A" y model because 
the progressive development of metabolic syndrome in yellow A" 3 
mice more accurately recapitulates the gradual onset of metabolic 
disease in humans with adult-onset obesity and diabetes. This is 
in contrast to genetic models of obesity and diabetes, which pro- 
duce early-onset metabolic disease," 5,18 or dietary manipulation 
where animals are forced onto a hypercaloric diet. 

A Western-style diet precipitates metabolic disease in 
obese-born offspring. In contrast to the overt defects in glucose 
metabolism seen in offspring of rodents fed a hypercaloric diet, 19 " 
21 mice born of obese yellow A vy la mothers showed no evidence 
of disturbed glucose homeostasis when compared with offspring 
of lean A vy ia mothers (Fig. 1; Fig. SI). However, offspring of 
obese mothers (ObC) were slightly but significantly heavier 
than offspring of lean mothers (LnC) before weaning, and males 
remained so afterwards (Fig. IB). We asked if this indicated a 
cryptic or latent predisposition to metabolic disease that might 
be unmasked by a dietary challenge. We exposed both groups to 
a Western-style high fat diet and found that while both groups 
(LnHF and ObHF) gained more weight on this diet, male ObHF 
mice became significantly heavier than LnHF by 12 weeks 
(Fig. IB). After only three weeks on the Western diet (at 6 weeks 
of age) male ObHF mice demonstrated glucose intolerance 
(Fig. 1C) and insulin resistance (Fig. ID); their female siblings 
displayed slight insulin resistance but no impaired response to 
glucose (Fig. 1C and D). Male (and to a lesser extent female) 
ObHF mice also exhibited a spike in blood glucose soon after 
intraperitoneal injection with insulin, indicating an impairment 
in their insulin response (Fig. ID). 

To further assess the metabolic phenotype induced by nat- 
ural-onset maternal diabetes and obesity, we sacrificed animals 
at 12 weeks of age for tissue analysis. We analyzed only male 
offspring as they demonstrated the most marked impairment in 
metabolic homeostasis. Liver lipidomic profiling revealed that 
total triacylglycerides (TAGs) and diacylglycerides (DAGs) were 
significantly higher in offspring of obese mothers, and levels were 
further elevated by the Western diet (Fig. 2A and B; Fig. S2). 
H&E stained sections of liver were scored using pathological cri- 
teria for nonalcoholic steatohepatitis in humans. 22 The Western 
diet induced steatosis, which was markedly more severe in ObHF 
mice (Fig. 2C and D) . 

We also assessed serum leptin, an adipokine that circulates at 
levels proportional to body fat. 23 Leptin was higher in ObC mice 
than controls (Fig. 2E) and was further elevated by the Western 
diet. It was also raised by the Western diet in offspring of lean 
mothers. Leptin correlated with body mass in all of the groups 
(Fig. 2F), suggesting that the observed increases in body weight 
are due to an increase in adiposity rather than an increase in lean 
mass. 

Taken together, these results indicate that the offspring of 
obese mothers have a latent defect in lipid and glucose metabo- 
lism that is unmasked by an energy-dense Western-style diet. 

Maternal obesity and diabetes alters hepatic gene expression 
in male offspring. The phenotypic defects in ObC and ObHF 
mice cannot have a genetic basis, since the mice are genetically 
identical to LnC and LnHF mice. In order to obtain a more 
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Figure 1. Offspring exposed to maternal obesity and diabe- 
tes have a latent predisposition to metabolic disease, which 
is revealed by a high-fat Western-style diet. (A) Experimental 
strategy. (B) Body weights of male and female offspring in 
each of the four groups in (A) (male: LnC n = 14, ObC n = 
13, LnHF n = 15, ObHF n = 9; female: LnC n = 19, ObC n = 13, 
LnHF n = 13, ObHF n = 18); inset, pre-weaning weights for 
Lean-born and Obese-born animals. (C) Blood glucose levels 
in male and female animals during a glucose tolerance test 
at 6 weeks of age (male: LnC n = 14, ObC n = 13, LnHF n = 
15, ObHF n = 9; female: LnC n = 20, ObC n = 6, LnHF n = 13, 
ObHF n = 11). (D) Blood glucose levels in animals during an 
insulin tolerance test at 7 weeks of age (male: LnC n = 14, 
ObC n = 17, LnHF n = 19, ObHF n = 19; female: LnC n = 20, 
ObC n = 13, LnHF n = 13, ObHF n = 18). Data are represented 
as mean ± SEM. Statistically significant differences are as 
indicated: *p < 0.05, **p < 0.01, ***p < 0.001 (LnC vs. ObC); 
+ p < 0.05, n p < 0.01, m p < 0.001 (LnHF vs. ObHF). 



detailed picture of the defects induced by maternal 
obesity and diabetes, we assessed hepatic gene expres- 
sion in 16 mice (4 from each group) using Affymetrix 
GeneChip Mouse Gene 1.0 ST arrays, which provide 
whole transcript coverage of 28,853 genes in the mouse 
genome. Using the statistical packages Messina 24 and 
LimmaGP (see Materials and Methods), we identified 
45 transcripts (42 unique genes) whose expression was 
significantly altered by maternal obesity (ObC vs. LnC; 
Fig. 3A) . Candidate gene expression changes identified 
by LimmaGP were, as expected, concordant across all 
animals within a group (because it uses a moderated 
t-statistic), whereas changes identified by Messina had 
larger effect sizes but more heterogeneity within groups 
(because it allows for intra-group variation); 24 Messina 
candidates were also more likely to be associated with 
a CpG island (Fig. 3A). qPCR was used to interrogate 
five candidate genes (chosen at random) associated with 
CpG islands, four of which exhibited altered expression 
in ObC relative to LnC as predicted by the array {Dbp, 
Gml29, Gnal4, Zbtbl6; Fig. S3). Validation included 
many animals that were not used on the array (17/22), 
and the same intragroup variation was observed in 
this larger cohort, particularly in the ObC group; this 
is consistent with a stochastic epigenetic response to 
maternal obesity. 

A similar pattern of expression changes was found 
in control mice exposed to the Western diet (LnHF vs. 
LnC); using the same array analysis parameters we find 
that 30 of the 45 transcripts (67%) identified as altered 
in ObC are also significantly altered in LnHF. This 
suggests that intrauterine and preweaning exposure to 
maternal obesity induces expression changes in a set of genes that 
are also responsive to overnutrition in later life. A much larger 
number of genes exhibited altered expression in the ObHF group 
compared with the other groups (Fig. S4); this is likely to reflect 
the obvious liver pathology in the ObHF mice (see Fig. 2C 
and D), which confounds any attempt to link gene expression 
changes to the latent predisposition in this group. 
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Gene ontology (GO) analysis of genes whose expression was 
altered in ObC as compared with LnC obtained a large set of 
ontologies, within which metabolic functions were present but 
not prominent (Fig. 3B). The two genes with the most strongly 
downregulated expression were ATPASE6 and CYTB, both 
encoded by the mitochondrial genome. We asked if this decrease 
reflected a lower mitochondrial copy number in hepatocytes of 
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Figure 2. Offspring exposed to maternal obesity and diabetes have defects in lipid metabolism which are exacerbated by a Western diet. Total liver 
(A) triacylglycerol and (B) diacylglycerol levels (LnC n = 4, ObC n = 6, LnHF n = 12, ObHF n = 8), and (C) steatosis scores (LnC n = 4, ObC n = 9, LnHF n 
= 10, ObHF n = 10) in livers of 12 week old male animals. Representative liver histology sections are shown in (D) (scale bar, 100 ^im). (E) Serum leptin 
levels and (F) relationship between serum leptin levels and body weight at 12 weeks (LnC n = 6, ObC n = 7, LnHF n = 12, ObHF n = 12). Box and whisker 
plots show the median, 25th and 75th percentile values with whiskers indicating the maximum and minimum; means are indicated by an X. Statisti- 
cally significant differences are indicated by asterisks: *p < 0.05, **p < 0.01, ***p < 0.001. Right-tailed Fisher's exact test, LnHF c.f. ObHF steatosis scores: 
p = 0.012. 



ObC mice, but found no difference in mitochondrial genome 
copy number between LnC and ObC (Fig. 3C). Among the most 
strongly upregulated genes are major urinary proteins (MUPs), 
which have been implicated in regulation of energy expenditure. 25 
Maternal obesity causes widespread epigenetic changes in 
male offspring. The characteristics of the mice exposed to mater- 
nal obesity and diabetes demonstrate that they are programmed in 
fetal life to respond to an energy-dense diet by developing adipos- 
ity, insulin resistance, glucose intolerance and hepatic steatosis. 
This programming may involve stable alterations in gene expres- 
sion states, mediated by epigenetic modifications. The molecular 
basis of epigenetic gene regulation is complex and incompletely 
characterized; however the role of cytosine methylation has been 
extensively investigated and is best understood in the suppression 
of transcription initiation. 26 Cytosine methylation in mammals 
is almost completely confined to CpG dinucleotides, which are 
concentrated in "CpG islands" that are transcriptional regulatory 
elements. 



We compared the genome wide patterns of CpG island meth- 
ylation in the livers of ObC and LnC mice (high fat-fed mice 
were not examined, as their overt metabolic disease would make 
it impossible to distinguish between methylation changes that 
were potentially causal of the latent phenotype, and those merely 
consequential to their disease). As previously described, 8 we 
made libraries enriched for the unmethylated fraction of DNA by 
sequential digestion with Hpall and McrBC, followed by ligation- 
mediated PCR. Libraries were hybridized to Agilent Mouse CpG 
Island 105K arrays covering approximately 16,000 CpG islands. 
We analyzed 8 male LnC mice and 8 male ObC mice; pooled 
libraries from 10 LnC controls acted as the reference sample for 
each array. 

We first used principal components analysis (PCA) to visu- 
alize the overall distribution of array data. PCA reduces many 
variables within a complex data set to a few artificial variables, 
which together account for most of the variance in the data. The 
first four principal components of our data accounted for over 
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Figure 3. The latent metabolic phenotype is associated with changes in hepatic gene expression. (A) Heat map of gene transcripts with significant 
expression changes in ObC vs. LnC offspring (n = 4 per group); note the similar patterning across all animals exposed to overnutrition (prenatally, 
postnatally or both). The columns on the right indicate genes that have CpG island promoters (green) or do not (white), and whether the gene was 
identified by LimmaGP (white), Messina (black) or both (gray). Average log 2 -fold expression changes are shown at the far right. (B) Bar graph show- 
ing -log p-values (with Benjamini-Hochberg correction for multiple testing) for gene ontologies significantly overrepresented in genes that are dif- 
ferentially expressed in ObC animals relative to LnC. A p-value cutoff of 0.05 is indicated by the red line. (C) Bar graph showing relative copy number of 
mitochondrial genes Atp6 and Nd6 (relative to the nuclear gene Ndufvl) in livers of LnC, ObC, LnHF and ObHF animals. Data are represented as mean ± 
SEM. Statistically significant differences are indicated by asterisks: *p < 0.05, **p < 0.01, ***p < 0.001 . 



half (54.4%) of the variability and are shown in two-dimen- 
sional score plots (Fig. 4A). The intermingling of ObC and LnC 
mice in principal components (PC) 1 and 2 indicates that most 
variation is common within these two groups. This spontane- 
ous inter-individual variation is reminiscent of the variation we 
have previously described in C57BL/6J mice. 8 However, distinct 
clustering and separation of ObC and LnC offspring is evident 
in PC 3 and 4. Taken together, the intermingling in PC 1 and 
2 and separate clustering in PC 3 and 4 indicate that intergroup 



variation resides within a larger context of the spontaneous varia- 
tion we and others 8,10 have previously described. A GO analy- 
sis of the 5% most methylation-variable loci in LnC and ObC, 
which are responsible for PC 1 and 2, produces a set of ontolo- 
gies related to developmental functions (Fig. 4B) that are nearly 
identical to the set identified in our previous study of intrinsic 
epigenetic variation in the C57BL/6J strain. 8 

We investigated the array probes responsible for the sepa- 
rate clustering of ObC and LnC in PC 3 and 4. Statistical 
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testing by LimmaGP revealed an enrichment of small ^-values 
(Fig. 4C), indicating that there are significant methylation dif- 
ferences between the two groups. Limma detected methyla- 
tion differences at 2350 probes corresponding to 1598 genes 
(at p < 0.01). Given the substantial intragroup variation, we also 
used Messina to identify stochastic changes in methylation that 
affected most but not all individuals in the ObC group. Most 
margin scores were very small (Fig. 4D); only 36 loci showed a 
margin of > 1.5 -fold change. The many more significant meth- 
ylation changes identified by LimmaGP (Fig. 4C) are therefore 
likely to be small in magnitude. Validation by bisulfite allelic 
sequencing showed that 6 of 9 loci examined exhibited small 
methylation differences (0.9—15.3%, Fig. S5) in ObC ani- 
mals in the direction predicted by the array. The verification 
rate (FDR -0.33) and the small magnitude of the changes are 
consistent with that of other studies of complex phenotypes 
using genome-wide strategies. 8,27 " 30 There was no significant 
overlap between the genes whose expression was changed in the 
ObC group and the genes whose methylation states differed. 
However, GO analysis of the genes identified by LimmaGP as 
differentially methylated (at p < 0.01) between the ObC and 
LnC groups again revealed significant enrichment in ontolo- 
gies associated with development (Fig. 4E); these are largely the 
same ontologies (8/10) as those in the intrinsically methylation- 
variable loci (Fig. 4B). Despite the prominent metabolic phe- 
notype in Ob mice, metabolic functions are not represented in 
these ontologies. This may suggest that the effects of maternal 
obesity are much broader than the metabolic aberrations uncov- 
ered by a Western diet. 

Discussion 

We have found that natural-onset maternal obesity and diabetes 
induces a latent defect in glucose and lipid metabolism in off- 
spring of A" y mice. This defect is unmasked by a Western-style 
diet, which causes mice to become insulin resistant and develop 
hepatic steatosis. Changes in gene expression induced by maternal 
obesity point to a broad set of functional ontologies, not limited 
to energy metabolism. Genome-wide assay of CpG island meth- 
ylation demonstrates that maternal obesity induces widespread 
methylation changes, which affect functional ontologies relating 
to development and gene expression. Thus the latent defect in 
energy metabolism is not prominently reflected in the expression 
and epigenomic changes induced by maternal obesity. 

The latent metabolic phenotype we observed in this study 
contrasts with those of other rodent models using high fat and 
sugar diets to mimic maternal overnutrition: in these cases glu- 
cose handling defects are apparent in offspring even when they 
are weaned onto a normal chow diet. 21,31 The syndrome of obese 
yellow A" y la mice more closely models the gradual onset and pro- 
gression of obesity and type 2 diabetes in humans. In this model, 
ala offspring exhibit only subtle defects in lipid metabolism and 
do not progress to a disease state if they are maintained on a stan- 
dard chow diet. Exposure to an energy-dense "Western" diet pro- 
duces overt effects, which result in mice exhibiting the equivalent 
of nonalcoholic steatohepatitis (NASH) in humans. If this effect 



occurs in humans, it would suggest that the children of obese 
and diabetic mothers are at risk of metabolic syndrome but can 
be protected by avoiding an energy-dense diet. 

Fetal programming by maternal nutrition has been proposed 
as a predictive-adaptive response to environmental conditions, 
whereby developing offspring adjust their physiology to suit the 
expected postnatal environment. While the underlying mecha- 
nisms are not clear, it is now generally supposed that epigenetic 
modification of metabolically relevant genes is involved. 32 But 
concrete examples are limited to a handful of genes, where the 
observed changes are small, 33 " 36 and very few exist in the context 
of maternal overnutrition. 37,38 Here we deliberately took an unbi- 
ased, genome-wide approach to ask whether epigenetic mecha- 
nisms are involved in the response to maternal obesity. While 
we found that maternal obesity induced methylation changes 
across the genome and expression changes at a variety of genes, 
metabolic genes were notable by their relative absence. Instead, 
the genes whose methylation is affected by maternal obesity 
cluster in developmental ontologies. This may indicate a greater 
susceptibility of developmental processes to environmental per- 
turbation, consistent with our previous finding that methylation 
patterns at genes in developmental ontologies are intrinsically 
more variable and susceptible to dietary supplements than other 
parts of the genome. 8 The systemic consequences of a pertur- 
bation in developmental processes could well underlie the 
metabolic phenotype and perhaps other cryptic phenotypes pro- 
grammed by maternal obesity and diabetes. Recent large-scale 
epidemiological studies on the prevalence of non-metabolic dis- 
orders, such as asthma and autism in children of obese moth- 
ers, support the idea that other defects may be programmed by 
maternal obesity. 39,40 

The methylation changes we observed were small in magni- 
tude, consistent with previous reports of epigenetic program- 
ming by us and others. 8,11,36,41,42 This raises the obvious question 
as to the relevance of the methylation changes to the phenotype. 
Certainly, we did not observe a large overlap between the genes 
with methylation changes and those with changed expression. 
This is reminiscent of the findings of Rando and colleagues, 
where widespread methylation changes in response to pater- 
nal high-fat feeding did not specify individual gene expression 
changes, but were nevertheless linked to a metabolic phenotype 
in offspring. 41 Notably, few metabolic genes were altered in this 
study, and those that were showed a magnitude of methylation 
change similar to that we observe here. 

The latency of the metabolic phenotype in our model permits 
the examination of methylation and gene expression changes that 
occur prior to the onset of disease, allowing us to avoid detect- 
ing changes that are consequential to disease. Consistent with 
their metabolic phenotype, the high-fat fed group exhibited gene 
expression changes related to metabolism; this is a common find- 
ing in fetal programming studies, and such changes are often 
assumed to be inborn. But in our control fed offspring, where 
disease is not yet present, we find that metabolic genes are prac- 
tically absent from comparisons between lean- and obese-born 
mice. The possible exceptions are the mitochondrially-encoded 
genes involved in oxidoreductase activity: mitochondrial activity 
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is known to play an important role in the pathogenesis of meta- 
bolic disease and particularly, in NASH. 43 We speculate that the 
changes we observe in expression of genes related to non-met- 
abolic functional ontologies may either play an indirect role in 
regulating metabolism, or indicate other cryptic defects in these 
mice, programmed by maternal obesity but possibly unrelated to 
energy metabolism. 

The predictive-adaptive view of fetal programming has been 
heavily influenced by the idea that the infants of undernour- 
ished mothers are adapted to an environment with scarce food 
resources. 1 An alternative possibility is that inappropriate nutri- 
tion (undernutrition, overnutrition or a modified supply of key 
nutrients) programs generalized defects that impair responses to 
environmental stresses or conditions. In this view, the relative 
paucity of metabolic genes in our lists of affected loci is not sur- 
prising. Since diet is always likely to be among the most promi- 
nent environmental conditions, defects related to the handling of 
nutrients may be the most readily unmasked; in our study this 
seems to particularly be the case with male offspring. Sexually 
dimorphic responses have been observed in a number of fetal pro- 
gramming studies (reviewed in ref. 44); while the effect of sex 
hormones is often invoked as the mediator of sexual dimorphism 
in the setting of metabolic programming, it may also be that 
sex-specific epigenetic profiles render either the male or female 
genomes more vulnerable to the effects of specific environmental 
insults, 45 either in utero or postnatally. 

The fetal response to maternal obesity is clearly maladaptive 
and the epigenomic response most likely represents some form 
of damage induced by a suboptimal intrauterine environment. 
However, the maladapted phenotype is not a fait accompli: 
When maintained on a healthy diet, the mice never develop met- 
abolic disease. Thus, although risk is programmed, overt disease 
is avoidable; if this translates to human populations, it highlights 
the need for multigenerational affirmative intervention to break 
the cycle of obesity. 

Materials and Methods 

Mouse breeding and diets. All animals were handled in accor- 
dance with good practice as defined by the National Health and 
Medical Research Council (Australia) Statement on Animal 
Experimentation and requirements of state government leg- 
islation. The study was approved by the Garvan/St. Vincent's 
Animal Ethics Committee (06/12 and 09/12). 

The A" y mice used in this study were descended from an iso- 
genic C57BL/6 colony at Oak Ridge National Laboratories and 
have been maintained at the Victor Chang Cardiac Research 
Institute since 2001. All parental mice were fed ad libitum on 
NIH-31 control diet (5% w/w fat, 13.5 MJ/kg). ala offspring 
were weaned onto either NIH-31 or SF00-219 high-fat diet 
designed to mimic a Western fast-food diet [equivalent to Harlan 
Teklad TD88137; 22% w/w fat (40% digestible energy), 0.15% 
w/w cholesterol, 19.4 MJ/kg]; A vy la littermates were culled at 
weaning. Feeds were manufactured by Specialty Feeds. 

Glucose and insulin tolerance tests. For metabolic testing, 
mice were fasted for 6 h and given an intraperitoneal injection 



of either 1.5 g glucose/kg mouse weight (for glucose tolerance 
test) or 0.75 units of insulin/kg mouse weight (for insulin toler- 
ance test). Blood glucose levels were then measured at 15 min 
(glucose) or 10 min (insulin) intervals using an Accu-Chek II 
Performa® glucose monitor (Roche Diagnostics GmbH). 

Lipid profiling. Lipidomic profiling was performed on liver 
samples as described previously. 46 

Liver histology. Haematoxylin and eosin stained sections 
were scored by an anatomical pathologist (M.E.B.) for nonal- 
coholic steatohepatitis according to standardized pathological 
criteria. 22 

Mitochondrial copy number. Mitochondrial copy number 
was assessed by qPCR quantification of mitochondrial genes 
Atp6 and Nd6 normalized to the nuclear gene Ndufvl. 47 Each 
reaction was done in triplicate with approximately 10 ng of 
DNA per reaction. See Table SI for primer sequences. 

Leptin assay. Serum was collected at 12 weeks, and leptin 
concentrations were measured using the Mouse Leptin ELISA 
(Millipore) following manufacturer's instructions. 

Statistical analysis for physiological data. All comparisons 
were made between age-matched, gender-matched, genetically 
identical mice. Statistical significance was determined by two- 
way ANOVA with two-tailed Student's t-tests for glucose and 
insulin tolerance tests, by two-tailed Student's t-tests for the 
other metabolic tests and by a right-tailed Fisher's exact test for 
comparison of steatosis scores. 

Gene expression. Total RNA was extracted from liver tis- 
sue and analyzed with GeneChip® Mouse Gene LOST arrays 
(Affymetrix). Labeling and hybridization were performed by 
the Ramaciotti Centre for Gene Function Analysis (University 
of New South Wales). Microarray data was preprocessed and 
normalized using robust multichip average (RMA), via the 
NormalizeAffymetrixST (version 2.0, available at pwbc. gar- 
van. unsw.edu. au/gp) GenePattern module (version 3.2.3 48 ). 
Differential gene expression analysis was performed using 
LimmaGP and Messina (see below). 

LimmaGP analysis. LimmaGP (version 19.3, available at 
pwbc. garvan.unsw.edu. au/gp; Cowley et al., manuscript in 
preparation) is a GenePattern module for identifying differ- 
ential gene expression using limma.^ Briefly, limma combines 
linear models to flexibly represent the experimental design of a 
microarray study in a mathematical framework, with an empiri- 
cal Bayes, moderated t-statistic, which has more power than the 
Student's t-test or ANOVA for analyzing microarray data. 49 We 
chose a cutoff of FDR < 0.05, with p < 1 x 10" 5 and a minimum 
fold change of 1.3, to capture genes with a consistent change 
between groups. 

Messina analysis. We used Messina 24 as a complimentary can- 
didate identification method. Messina takes a machine learning 
approach and allows for intra-group heterogeneity by permitting 
a user-specified, low level of misclassification. This allows detec- 
tion of candidate genes with inconsistent changes in expression, 
so is particularly suited to studies where changes may be sto- 
chastic. Messina represents the magnitude of differential expres- 
sion by a margin score (representing the difference in expression 
between the closest included members of the two groups). We 
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Figure 4. Livers of obese-born males exhibit widespread changes in cytosine methylation. (A) Two-dimensional plots showing principal component 
(PC) scores of LnC and ObC animals (n = 8 per group). Three plots are shown: scores for PC 1 vs. PC2 (left), PC2 vs. PC 3 (middle) and PC 3 vs. PC 4 
(right). The amount of variability accounted for by each PC is indicated on the axes. The wireframe ellipses surrounding PC scores for each group are 
indicative of the overall standard deviation within each group. (B) Bar graph showing -log p-values (with Benjamini-Hochberg correction for multiple 
testing) for gene ontologies significantly overrepresented in methylation variable genes in LnC and ObC animals. A p-value cutoff of 0.05 is indicated 
by the red dashed line. (C) Histogram showing the distribution of LimmaGP p-values across all probes on the methylation array. A p-value cutoff of 0.01 
is indicated by the red dashed line. (D) Histogram showing the distribution of Messina margin scores from the methylation array. Inset shows enlarged 
view of the histogram tails. (E) Bar graph showing -log p-values (with Benjamini-Hochberg correction for multiple testing) for gene ontologies identi- 
fied as being significantly overrepresented in genes differentially methylated in ObC vs. LnC (as detected by LimmaGP, p < 0.01). A p-value cutoff of 
0.05 is indicated by the red dashed line. 
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set a maximum misclassification of one animal out of four in 
both control and test groups (i.e., specificity and sensitivity 
equals 0.75), and applied a margin of > 0.585 (equivalent to 
> 1.5-fold change). 

Microarray validation. One microgam total RNA was reverse 
transcribed and assayed by qPCR. Each reaction was done in 
triplicate with approximately 10 ng of cDNA per reaction. 
Samples were normalized to the geometric mean of four house- 
keeping genes {Actb, Hprt, Gapdh and Canx); see Table SI for 
primer sequences. 

Methylation analysis. DNA was extracted, enriched for 
unmethylated sequence by Hpall and McrBC digest, and 
hybridized to two-color Agilent Mouse 105K CpG island arrays 
as previously described. 8 Each sample was hybridized against a 
reference pool made by combining DNA libraries from 10 LnC 
offspring. Data was analyzed with LimmaGP and Messina algo- 
rithms as described above. Principal components analysis was 
performed in Partek Genomics Suite (v6.5). 

Bisulfite sequencing validation. Allelic methylation patterns at 
nine candidate loci were assessed by bisulfite allelic sequencing. 50 
Two micrograms DNA was treated with sodium bisulfite using 
the Epitect kit (QIAGEN GmbH), PCR amplified, cloned into 
E. coli and sequenced. See Table SI for primer sequences. 



Gene ontology. Gene ontology analysis was performed on 
gene expression and methylation candidates using the standard 
workflow in Ingenuity Pathways Analysis (Ingenuity Systems). 
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